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Abstract. Grids of stellar atmosphere models and associated synthetic spectra are numeri- 
cal products which have a large impact in astronomy due to their ubiquitous application in 
the interpretation of radiation from individual stars and stellar populations. 3D model at- 
mospheres are now on the verge of becoming generally available for a wide range of stellar 
atmospheric parameters. We report on efforts to develop a grid of 3D model atmospheres 
for late-type stars within the CIFIST Team at Paris Observatory. The substantial demands 
in computational and human labor for the model production and post-processing render this 
apparently mundane task a challenging logistic exercise. At the moment the CIFIST grid 
comprises 77 3D model atmospheres with emphasis on dwarfs of solar and sub-solar metal- 
licities. While the model production is still ongoing, first applications are already worked 
upon by the CIFIST Team and collaborators. 
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1. The simulation code 

The 3D simulations were performed 
with the radiation-hydrodynamics code 
CQ 5 BOLD dFreytag. Steffen. & Dorchl l2002t 
IWedemever et al.l 2004I) . The code solves the 
time-dependent equations of compressible 
hydrodynamics coupled to radiative transfer 
in a constant gravity field in a Cartesian 
computational domain which is representative 
of a volume located at the stellar surface. The 
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equation of state takes into consideration the 
ionization of hydrogen and helium, as well as 
the formation of H2 molecules according to 
Saha-Boltzmann statistics. Relevant thermo- 
dynamic quantities - in particular gas pressure 
and temperature - are tabulated as a func- 
tion of gas density and internal energy. The 
multi-group opacities used by C0 5 BOLD are 
based on monochromatic opacities stemming 
from the MARCS stella r atmosphere package 
(IGustafsson et al.l I2008I) provided as function 
of gas pressure and temperature with high 
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wavelength resolution. For the calculation 
of the opa cities solar element al abund ances 
according iGrevesse & Sauvall (119981) are 
assumed with the exception of CNO for 
w hich valu e s clos e to the recommendation 
of lAsplundl d2005l) are adopted (specifically, 
A(C)=8.41, A(N)=7.80, A(0)=8.67). 



2. Model set-up 

We typically used a number of 140 x 140 x 150 
points for the hydrodynamical grid. The 
wavelength dependence of the radiation 
field was represented by 5 multi-group 
bins in the case of solar metallicity, and 6 
bins at sub-solar metalliciti es follow i ng th e 
procedu res laid o ut by |Nordlundl ( 1982); 
iLudwid dl992l); ILudwig. Jordan & S teffen 
d 19941) : IVogler. Bruls. & Schiisslerl d2004l) . 
For test purposes we calculated a few 
models using 12 bins. The sorting into 
wavelength groups was done applying thresh- 
olds in logarithmic Rosseland optical depth 
{+00,0.0,-1.5, -3.0,-4.5,-oo} for the 5- 
bin, and {+oo, 0.1, 0.0, -1.0, -2.0, -3.0, -oo} 
for the 6-bin schemes. In the case of 12 
bins we used as thresholds {+oo, 0.15, 0.0, 
-0.75, -1.5, -2.25, -3.0, -3.75, -4.5, -oo}; 
in addition, the first three continuum-like 
bins were split each into 2 bins according to 
wavelength at 550, 600, and 650 nm. In all 
but one bin a switching between Rosseland 
and Planck averages was performed at a band- 
averaged Rossland optical depth of 0.35; in 
the bin gathering the largest line opacities the 
Rosseland mean opacity was used throughout. 
The decisions about number of bins, and 
sorting thresholds are motivated by comparing 
radiative fluxes and heating rates obtained by 
the binned opacities in comparison to the case 
of high wavelength resolution. 

We constructed 3D starting models by scal- 
ing existing 3D models according to ID stan- 
dard stellar structure models. The ID models 
provided scaling factors for thermodynamic, 
kinematic, and geometrical properties to mod- 
ified r e ff, logg, and metallicity. Points of the 
same optical depth and same pressure rela- 
tive to the surface pressure were considered 
associated between the original and scaled 



model, in the optically thin and thick regions, 
respectively. Empirically, the scaling proce- 
dure works particularly reliably if the scal- 
ing is performed along lines of constant en- 
tropy jump in the r B ff-log g-pla ne (see, e.g., 
ILudwig. Frevtag.& Steffenlll999l) . 

At given T e ff and log g we kept the compu- 
tational domain the same for each metallicity. 
Since the models were intended to serve as at- 
mosphere models for spectral synthesis calcu- 
lations, a minimum extent of the computational 
domain to log Tr oss = -6 was always required. 
The overall philosophy was to keep the model 
set-up - in a relative sense - solar-like. We re- 
quired a minimum time interval over which the 
model could be considered as relaxed which 
corresponded to an equivalent of one hour in 
a solar model. 



3. Model calculations and monitoring 

The majority of the simulation runs was con- 
ducted on a dedicated Linux cluster of 14 
double-processor, double-core nodes running 
one model per node making use of OpenMP 
parallelization. A main-sequence run is com- 
pleted in 1-2 months (wall clock time). More 
CPU demanding giant models were calculated 
on an IBM SP5 at the Italian CINECA super- 
computing center, and in addition on machines 
at the Institute of Theoretical Physics and 
Astronomy of Vilnius University (Lithuania). 
Running up to 20 hydrodynamical models in 
parallel made it necessary to facilitate the mon- 
itoring of simulation runs. For this purpose 
we developed an automatic procedure produc- 
ing diagnostic plots in HTML-format depicting 
the evolution of the thermal, flux, granulation 
pattern, and oscillatory properties of a model. 
Figure Q] shows examples of granulation pat- 
terns at different metallicity. 

4. Snapshot selection and auxiliary 
data. 

After completion of a model run a sub-set 
of snapshots is selected for spectral synthesis 
purposes. The selection is guided by the re- 
quirement that the statistics of the sub-sample 
should closely resemble the statistics of the 
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Fig. 1. Typical snapshots of the emergent bolometric intensity (false color representation) from two runs 
at r eff =5000K, logg=2.5 and metallicties [M/H]=0 and -2. The metallicity, the relative RMS intensity 
contrast of the image and a scale of 10 times the pressure scale height at r = 1 are given in each frame. 



whole run. In particular, the statistics of fluc- 
tuation in velocity and emergent flux should 
be preserved. The selection is done for reduc- 
ing the amount of necessary calculations since 
3D spectral synthesis calculations are com- 
putationally demanding. The selected data is 
augmented by basic information about mean 
model properties, and ID standard atmosphere 
models which exactly share the atmospheric 
parameters with the 3D run. 



5. Grid completeness and 
computational costs 

At the moment 77 non-solar 3D models have 
been completed (see Fig. [5}. On top of this 
number we calculated a series of solar mod- 
els and test models with more detailed (12 
bins) radiative transfer. Due to a combina- 
tion of specific needs for projects and com- 
putational costs most models were produced 
for dwarfs and sub-giants. Figure [3] illus- 
trates the approximate computational costs - 
i.e., the demand in CPU time - of a model 
run; the symbol diameter is proportional to 
the quantity (Aj/A?cfl) -0 ' 5 where Af is the 
typical timestep of a model run, and A?cfl 
the Courant-Friedrichs-Levy time. Metallicity 



generally plays no dominant role for the com- 
putational cost. The computational cost rises 
significantly towards the red giant branch, and 
towards higher effective temperatures. This ba- 
sically reflects the decrease of the radiative 
time scale relative to A?cfl leading to a reduc- 
tion of the overall time step. 



6. Thoughts for a next generation 3D 
grid 

At present the realism of 3D model atmo- 
spheres is limited by the approximations nec- 
essary to keep the radiative transfer tractable. 
Consequently, in next generation grids much 
effort will be put in refining the radiative trans- 
fer. This will be accompanied by a - likely 
modest - increase in spatial resolution. A per- 
haps less obvious property should be a higher 
degree of homogeneity among the models. A 
global strategy for designing the computational 
grid is desirable which takes into consideration 
the physical structure and constraints given by 
the numerics. This also relates to the way opac- 
ities are grouped. A global strategy ensures 
an increased differential accuracy of proper- 
ties among models (e.g., granulation contrast, 
granular scale, equivalent mixing-length pa- 
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Fig. 2. Status of the model production for the CIFIST grid. Symbols mark the location of a model in the 
r c ff-log g-plane. Green color indicates completed model runs, red ongoing calculations. Each square is split 
into four sub-squares indicating solar, 1/10, 1/100, and 1/1000 of solar metallicity. The solar position is 
indicated by the round symbol. 



rameter, equivalent width of synthetic lines, 
etc.). This process will benefit from the already 
existing 3D models. If grid completeness is im- 
portant the computational cost of giants neces- 
sitates to include them in the production early. 
Last but not least dedicated efforts are neces- 
sary to make the computational results acces- 
sible to the broader astronomical community. 
Extensive tables of synthetic colors and spec- 
tra are the obvious data products which need to 
be delivered. 
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Fig. 3. Approximate computational model cost. 
The surface area of a circles depicts the approximate 
compute time of a model run. For details see text. 
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